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Abstract Alfvenic waves are thought to play an im- 
portant role in coronal heating and solar wind accelera- 
tion. Recent observations by Hinode /SOT showed that 
the spicules mostly exhibit upward propagating high 
frequency waves. Here we investigate the dissipation of 
such waves due to phase mixing in stratified environ- 
ment of solar spicules. Since they are highly dynamic 
structures with speeds at about significant fractions of 
the Alfven phase speed, we take into account the effects 
of steady flows. Our numerical simulations show that 
in the presence of stratification due to gravity, damp- 
ing takes place in space than in time. The exponential 
damping low, exp(-At 3 ), is valid under spicule condi- 
tions, however the calculated damping time is much 
longer than the reported spicule lifetimes from obser- 
vations. 
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1 Introduction 



The coronal heating (|Edlenlll943l ) mechanism is one of 
the major unsolved problems in solar physics. Since the 
energy flux carried by acoustic waves is too small, the 
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possibility of heating by MHD waves has been investi- 
gated intensively as the magnetic structure of the solar 



coron a can play an important role here (jHood et al 



1997|) ■ The propagation of Alfven waves is one of the 



candidate mechanisms that can carry energy to large 
distances from the surface, and heat the solar corona. 
However, a heating theory based on the Alfven waves 
faces a couple of difficulties: Firstly, the waves have 
to transport enough energy flux, and secondly, they 
have to dissipate efficiently in order t o deposit the right 



amou nt of energy at the right place ()De Moortel et al 
1999). The Alfven waves may reach the corona even 



in the absence of highly stratified atmosphere but 
with lesser propagation speed. The damping length 
of Alfven modes is defined b y various dissipative pro- 
cesses such as p hase mixing ( Hevvaerts fc Priestlll983t 
Browning 1991 ). resonance abs orption (llonson 19781) . 



and nonlinear mode conversion ( Hollweeill982l ). Phase 



mixing is a mechanism for dissipating Alfven waves , 
which was first proposed by iHevvaerts fe Priest! (j 19831 ) . 
When Alfven waves propagating in an inhomogeneous 
medium, on each magnetic field line, a wave propagates 
with its own local Alfven speed. After a certain distance 
or after enough time, these neighboring perturbations 
will be out of phase. This ultimately results in a strong 
enhancement of the dissipation of Alfven waves energy 



via both viscosity and resisti vity. iHood et al.1 (|2002l ) 
and IHevvaerts fc Priestl (| 19831 ) analytically showed that 
in both the strong phase mixing limit and the weak 
damping approximation, the ampl itude of Alfven waves 
decay s with time as exp(-t 3 ). iKarami fc Ebrahimi 
(2009) calculated numerically the damping times of os- 
cillations in the presence of viscosity and resistivity in 
coronal loops. They concluded that the above expo- 
nential damping law in time is valid for the Lundquist 



numb ers higher than 10 7 . iDe Moortel et al 



d 19991 



2000l) studied the effect of stratification due to grav 



ity on phase mixing, and found that the wavelengths 



2 



lengthen when Alfven waves propagate through a strat- 
ified plasma. They concluded that a vertical stratifica- 
tion of density makes phase mixing by ohmic heating 
less important in coronal heating problem. They also 
found that in a stratified atmosphere, the heat will be 
deposited higher up than in an unstratificd atmosphere, 
and that the viscous heating will be the dominant com- 
ponent in the heating processes at lower heights. How- 
ever, in coronal conditions the effect of stratification on 
efficiency of phase mixing would still be large as the 
height at which most heat would be deposited through 
ohmic dissipation is increased considerably by stratifi- 
cation. Moreover, depending on the value of pressure 
scale height, phase mixing can either be more or less 
efficient than in the uniform case. 

Spicules ha ve long been investigated a s a coronal 
heating agent (jrlollwe d [l983 lAthavlll983) . They are 
grass-like spiky features seen in chromospheric spectra l 
lines at the solar limb (jZaaarashvili &: Erdelvil 120091 ) . 



These spiky dynamic jets are propelled upwards at 
speeds of about 20 — 25 km s _1 from photosphere into 
the magnetized low atmosphere of the sun. Their di- 
ameter varies from spicule to spicule having the values 
from 400 km to 1500 km. The mean length of classi- 
cal spicules varies from 5000 km to 9000 km, and the 
typical life time of them is 5 — 15 min. The typical 
electron density at heights where the spicules are ob- 
served is 3.5 x 10 10 — 2 x 10 11 cm -3 , and th eir tem- 
perat ure s are estimated ~ 5000 — 8000 K (IBeckers 



19681) . Kukhianidze et all ( 20061 ); Zaaarashvili et al 



(2007|) by analyzing the height series of Ha spectra 
in solar limb spicules observed their transverse oscil- 
lations. The period of them estimated 20 — 55 and 
75 — 110 s. They concluded that these oscillations 
can be caused by propa gating kink waves in spicules. 



De Pontieu et al.l (|2007l) based on Hinode observations 
concluded that the most expected periods of transverse 
oscillations lay between 100 — 500 s, which interpreted 



as sign atures of Alfven waves. Okamoto fc De Pontieul 
( 20111 ) used Hinode /SOT observations of spicules, and 



concluded that upward propagating, downward prop- 
agating and standing waves occurred at the rates of 
about 60%, 20% and 20%, respectively. Furthermore, 
they found that upward propagating waves dominate 
at lower latitudes, and the medians of amplitude, pe- 
riod, and velocity amplitude are 55 km, 45 s, and 7.4 



km s _1 , respectively. More recently lEbadi et al.l (|2012l) 
made time-slice images of spicules, which were observed 
by Hinode/SOT. They concluded that the energy flux 
stored in spicule axis oscillations is of the order of coro- 
nal energy loss in quiet sun. These results motivated 
us to study the phase mixing of upward propagating 
Alfven waves in a stratified atmosphere. To do so, the 



section 2 gives the basic equations and the theoretical 
model. In section 3 the numerical results are presented 
and discussed, and a brief summary is followed in sec- 
tion 4. 



2 Theoretical modeling 

In the present work we keep the effects of stratification 
due to gravity in 2D x-z plane. The phase mixing and 
dissipation of propagating Alfven waves in a region with 
non-uniform Alfven velocity is studied. MHD equations 
governing the plasma dynamics are as follows: 



dp 



dv 



v ■ (H = o, 



(1) 



pg+ -(V x B) x B + pv\7 2 v, 

(2) 



- = Vx(vxB) 



pRT 



t?V 2 B, (3) 



(4) 



V • B = 0, 



(5) 



where v and rj are constant viscosity and resistivity co- 
efficients, and other quantities have the usual mean- 
ing. In particular, typical values for rj in the solar chro- 
mosphere and corona are 8 x 10 8 T~ 3 / 2 and 10 9 T~ 3 / 2 
m 2 s -1 , respectively. The value of pv for a fully ionize d 



( Priestlll982h 



H plasma is 2.2 x Kr 17 T 5 / 2 kg 
We assume that the spicules are highly dynamic with 
speeds that are significant fractions of the Alfven speed. 
The perturbations are assumed independent of y, with 
a polarization in y direction, i.e.: 



v = v k + v y (x, z,t)j 
B = Bgk + by(x, z,t)j . 



(0) 



Therefore, the pressure gradient is balanced by the 
gravity force, which is assumed to be g=g k via this 
equation: 



-Vp + Pog = 0, 



(7) 



3 



and the pressure in an equilibrium state is: 
Pa =p (x)e- z/H . 

The density profile is in the form of: 

po = po(x)e~ z/H , 
with 

RT 



H 



(8) 



(9) 



(10) 



where H is the pressure scale height. The linearized 
dimcnsionless MHD equations with these assumptions 
are: 



Fig. 1 The plot of initial wave packet with d=0.3a is 
presented. 



8Vy 8Vy ~. .dby 

and 



(11) 



v y (x = 0, z, t) — v y (x — 2, z, t) = 
b y (x = 0, z, t) = b y (x = 2, z, t) = 0, 



(15) 



db 
~dt 



az 



8Vy 

dz 



(12) 



where the velocities, the magnetic field, time and space 
coordinates are normalized to Vao = Bo/ 'yfixpo (with 
po as the plasma density at z = 0), Bq, t& (the pe- 
riod of Alfvcn waves), a (spicule radios), respectively. 
Also the resistivity and viscosity coefficients are nor- 
malized to a 2 /tA- The second terms in the left hand 
side of eqs.Qj] and fT2l present the effect of steady flows. 
Va{x, z) is the Alfven velocity, which for a phase mixed 
an d stratified atmosphere due to gravity is assumed t o 
be (|De Moortel et al.lll999HKarami fc Ebrahimill2009h : 



V A (x,z) = V A0 e z/2H [2 + t&nh[a(x- 1)]], 



(13) 



where parameter a controls the size of inhomogcneity 
across the magnetic field. 

The set of eqs.QTJ and [T2l should be solved under these 
initial and boundary conditions: 



v y (x, z, t = 0) 
b y (x, z, t = 0) 



Vao exp 



1/Z - ^2 
~2 { d 



-r 



Asm(nx) sin(7rz), 



sm(kz)e z/4H 
(14) 



where d is the width of the initial packet and A = 
10~ 7 . Figured] is the plot of initial wave packet given 
by equation 1141 for d = 0.3a (a is the spicule radius). 
The parameter k is chosen in such a way to have upward 
propagating Alfvcn wave. 



dvy 

at 



o. 



(16) 



3 Numerical results and discussion 

To solve the coupled eqs. [TTJ and [T^] numerically, the 
finite difference and the Fourth-order Runge-Kutta 
methods are used to take the space and time deriva- 
tives, respectively. We set the number of mesh grid 
points as 256 x 256. In addition, the time step is cho- 
sen as 0.002 (time is normalized to the Alfven time, 
t A ), and the system length in the x and z dimen- 
sions (simulation box sizes) are set to be 2000 km and 
6000 km. The parameters in spicule env i ronment are 
as fol lows ( Murawski fc Zaqarashvili 2010l : Ebadi et al. 
l2012h : a (spicule radios)=500 km, d = 0.3a = 150fcm 
(the width of gaussian paket), L=6000 km (Spicule 
length), v = 25km/s, B Q = 10G, n e = 10 n cm~ 3 , 
T = 8000^, g = 272ms- 2 , R = 8300m 2 S - 1 k- 1 
(universal gas constant), Vao = 40/cto/s, k = An/3 
(dimcnsionless wavenumber normalized to a), v = 
10 3 m 2 s _1 , T) = W 3 m 2 s- 1 , pi = 0.6, H = 500km, 
\a= 37.5 s (the period of Al fven waves), and a = 2 
( Okamoto fc De Pontieu1l201ll ). 



Figure [2] shows the perturbed velocity variations 
with respect to time in x = 1000 km, z = 1000 km; 
x = 1000 km, z = 3000 km; and x = 1000 km, z = 5000 
km respectively. We presented the perturbed magnetic 
field variations obtained from our numerical analysis in 
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Figure [3] for x = 1000 km, z = 1000 km; x = 1000 
km, z = 3000 km; and x = 1000 km, z = 5000 km 
respectively. In these plots the perturbed velocity and 
magnetic field is normalized to Vao and Bq respectively. 
In each set of plots it is appeared that both the per- 
turbed velocity and magnetic field are damped in the 
developed stage of phase mixing. 
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Fig. 3 The same as in Figure [2] but for the perturbed 
magnetic field. 
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Fig. 2 The perturbed velocity variations with respect to 
time in x = 1000 km, z = 1000 km; x = 1000 km, z = 3000 
km; and x = 1000 km, z = 5000 km respectively from top 
to bottom are showed. 
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Further, at the first height (1000 km), total ampli- 
tude of both velocity and magnetic field oscillations 
have values near to the initial ones. As height in- 
creases, the perturbed velocity amplitude does increase 
in contrast to the behavior of perturbed magnetic field. 
Nonetheless, exponentially damping behavior is obvi- 
ous in both cases. This means that with an increase in 
height, amplitude of velocity oscillations is expanded 
due to significant decrease in density, which acts as in- 
ertia against oscillations. Similar results are observed 
by time-distance analysis of Solar spicule oscillations 
( Ebadi et al]|2012l ) . It is worth to note that the density 
stratification influence on the magnetic field is negligi- 
ble, which is in agreement wit h Solar Optical Tel escope 
observations of Solar spicules (jVerth et al.ll201ll) . 

Figures 0] [5] show the 3-D plots of the perturbed 
velocity and magnetic field with respect to x, z for 
t = 30tA, t = 60tA, and t = 80^. They show that in 
the presence of stratification due to gravity, the damp- 
ing takes place in space than in time as an important 
point of these graphs. It should be emphasized that 
the damping time scale of the velocity field pattern 
is longer than the corresponding magnetic field pat- 
tern due to the initial conditions. In other words, in 
spite of standing waves, propagating waves are stable 
and dissipate after some periods due to phase mixing 
( Hevvaerts fe Priest|[l983h . 

For calculating the damping time it is suitable to 
calculate the total energy (kinetic energy plus magnetic 
energy) per unit of length in y direction as: 



E t ot{t) 



16tt 



dx 



dz[p(x, z)vy(x, z) + by(x, z)]. 

(17) 



In Figure |5] we plot the normalized total energy for 
three initial wave packet widths, i.e. d = 0.1a, d = 0.3a, 
and d = 0.8a. 

Since the treatments of different profiles of total en- 
ergy are similar and the damping times are very close 
to each other, hence, we continue our calculations with 
d = 0.3a, which is more logical compared with tube 
radius. In Figure [7] the kinetic energy, magnetic en- 
ergy, and total energy normalized to the initial total 
energy are presented respectively from top to bottom. 
The damping time calculated from total energy pro- 
file is 1050 s. Actually, spicules have short lifetimes, 
and are transient phenomena. We claim that in such 
circumstances, phase mixing can occur in space rather 
than in time. This is in agreement with previous works 
i De Moortel et al.lll999l ). We can definitely consider it 
as an uncertain rule of thumb. 







Fig. 4 (Color online) The perturbed velocity in x — z space 
is presented. The panels from top to bottom correspond to 
t = 30tA, t = G0t a, and t = 80£a, respectively. 



The total energy profile is best fitted to the ex- 
ponential damping function in time, exp(-At s ), with 
A = 10~ 9 , B = 2.85. Since the spicules are structures 
with low resistivit y and viscosity coefficients , so this is 
in agreement with iHevvaerts fc Priestl (|1983h work. 

So far the presented results in this paper are taken 
with considering a steady flow (vo = 25km /s). Actu- 
ally, it is of interest to investigate the differences in our 
numerical analysis outputs if the flow effect is omitted. 
To do that, we repeated the numerical analysis with 
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Fig. 5 (Color online) The same as in Figure[4]but for the 
perturbed magnetic field. 

zero-steady flow. Interestingly, the relevant results did 
not show significant difference in respect to the case of 
non-zero steady flow. A short theoretical scaling may 
help to clarify this result; comparing two dimensionless 
terms in cq. 1111 the convective term, which is associ- 
ated with the steady flow, vo^jr, and the Lorentz term, 
V2(x,z)%*, yields: 
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Fig. 6 The normalized total energy for three different ini- 
tial wave packet widths is presented. It is normalized to its 
value at t — 0. 



~ e z/H by/v y , (18) 

Here b y , v y are the average values of b y , v y normal- 
ized to Bq and Vao- Also, vq is normalized to Vao as- 
suming that both have the same order of magnitudes. 
Moreover, Va ~ Vao& z ^ 2H ^ has been taken into ac- 
count. Now with the large values of z, and assuming 
that b y and v y have a same order of magnitudes, then 

e _ k» ^> i. Meaning that the effect of the Lorentz 
term which is associated with the stratification effect 
is more important than the dominated convective term 
bringing into account the steady flow effects. Thus as 
long as z/H S> 1, one might ignore the steady flow 
effects here. Similar scaling can be done in eq. Q21 to 
show that in the presence of stratification the steady 
flow effects arc ignorablc. 

4 conclusion 

In our simple model, we assume that spicules are small 
scale structures (relative to coronal loops and other 
mega structures) with a uniform magnetic field, and 
a uniform temperature in all heights, and the transi- 
tion region between chromosphere and corona has been 
neglected. Density change along spicule axis is consid- 
erable, and stratification due to gravity is significant. 
As a result, the medium is dense in its lower heights, 
but it becomes rare and rare as height increases. Also, 
spicules have short lifetimes, and are transient phenom- 
ena. We claim that in such circumstances, phase mix- 
ing can occur in space rather than in time. This is in 
agreement with previous studies. In ordinary heights of 
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Fig. 7 Time variations of normalized kinetic energy, mag- 
netic energy, and total energy for d = 0.3a are presented 
from top to bottom respectively. 



a spicule, there is enough space to amplitude of a wave 
to damp, and to transfer its energy to the medium. Fur- 
thermore, we observed that the amplitude of perturbed 
velocity field, magnetic field and total energy decreases 
with height exponentially. Our numerical analysis show 
that the main phase of evolution dynamics occurs only 
in the first one third of the spicule height. Furthermore, 
we repeated the numerical analysis with zero-steady 
flow. Interestingly, the relevant results did not show 
significant difference in respect to the case of non-zero 
steady flow. 
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